Quadratic Solitary Waves in a Counterpropagating Quasi-Phase-Matched 

Configuration 
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We demonstrate the possibility of self-trapping of optical beams by use of quasi phase matching in a 
counterpropagating configuration in quadratic media. We also show the predominant stability of these 
spatial self-guided beams and estimate the power level required for their experimental observation. 
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Self-guided optical beams (or spatial solitons) have at- 
tracted significant research interest because they offer the 
possibility of all-optical switching and controlling light 
by light (see, e.g., (j). During the past few years it 
has been realised that quadratic nonlinearity is particu- 
larly attractive for potential practical realizations of all- 
optical switching, in that it not only supports stable soli- 
tons both in planar waveguides and bulk media but also 
provids an ultra-fast electronic nonlinear response (see, 
e.g., |2|). The advantages of quadratic nonlinear ma- 
terials are hampered in part by difficulties in obtaining 
close phase- velocity matching between interacting waves. 
One of the most effective ways to achieve such matching 
is the use of the so-called quasi-phase-matching (QPM) 
technique, in which large wave- vector mismatch between 
interacting waves is compensated for by periodic alterna- 
tion of the sign of effective x^ 2 -* coefficient. This technique 
has been known since 1962, Ref. 0], but only in the past 
decade has technological progress put the QPM technique 
in the front line of modern nonlinear optics 0] . In spite 
of all the theoretical and experimental progress achieved 
in the field of quadratic solitons, only solitons formed by 
waves with the same direction of propagation have been 
analysed so far (conventional co-propagating configura- 
tion; see, e.g., @f|). In a few works the advantages and 
implications of QPM technique for this type of solitons 
have been analysed specifically |6j. However, a paramet- 
ric interaction between counterpropagating quasi-phase- 
matched waves in quadratic [x ] media is also possible. 
Corresponding analysis made for for non-soliton (plane 
wave interaction) case R] revealed certain advantages of a 
counterpropagating interaction system over conventional 
copropagating strategies. Moreover, recently a very simi- 
lar so-called backward QPM configuration has been inves- 
tigated experimentally (see ||). In this work we investi- 
gate the QPM counterpropagating scheme, searching for 
solitary waves and investigating their stability. 

We consider the interaction between four optical waves 
in a slab waveguide with an appropriate nonlinear grat- 
ing. Two forward-propagating waves, the fundamental 
at frequency u) and with wave number k^ and the sec- 
ond harmonic (2a;, k 2ul ), are coupled with two backward- 
propagating ones, the fundamental (a;, — k u ) and the sec- 
ond harmonic (2a;, — k 2w ) [see Fig. 61 (a)]. Spatial mod- 
ulation of nonlinear susceptibility along a crystal can 
be described in terms of square-wave function d{z) [see 
Fig. [H (b)]. In this case the only nonzero matrix el- 



ements of the Fourier transform of d(z) are given by: 
di = -2i/(nl), where I = 2m - 1, m = ±1,±2,±3... 
Following the method developed in Ref. M we can derive 
the corresponding normalised system of equations which 
has the following dimensionless form: 
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E^(x, z), E2 LU (x, z) are the envelopes of the fundamen- 
tal wave and its second harmonic, respectively, sign "+" 
("-") corresponds to the forward (backward) propagat- 
ing wave; x is the transverse coordinate normalised on the 
width of the beam ro; z is the propagation distance which 
is normalised on the diffraction length Zd = k^r 2 ,; pa- 
rameters a + and a_ are nonlinear induced propagation 
constant shifts of the fundamental waves. Other two sys- 
tem parameters, a = k 2u /k u — 2 and S = (2ttI/L— k 2uj )zd 
(where L is the period of the nonlinear grating and / is 
the order of QPM) , are defined by the particular experi- 
mental setup. Note, that in contrast to Ref. jq| we have 
omitted all effective cubic terms in Eqs. (pj). This is well 
justified for lower order QPM (for I < 15) because in this 
case the cubic terms would become noticeable only when 
the light intensity exceeds the damage threshold for the 
typical nonlinear crystal. 

The system (|l|) has the following family of power-like 
integrals of motion: 



Q( Pl ,p 2 )=f:™{2 Pl \E+\ 2 + 2p 2 \E- 
+a( Pl -p 2 )(\E+J 2 -\E^\ 2 )}dx, 



(2) 



where pi, p 2 are any real numbers. Using the asymptotic 
expansion technique and the method based on integrals 
of motion (see Ref. B) one can demonstrate that the sta- 
bility threshold for the fundamental family of stationary 
localised solutions of the system (fil) is given as: 



d(Qi,Q 2 ) = 0Q 1 dQ 2 dQ 1 dQ 2 
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where Q\, Q2 are any two linearly independent invari- 
ants from the family (|2j), calculated for the fundamental 
stationary solitons, e.g. Qi = Q(1,0) and Qi = Q(0, 1). 
More elaborate analysis shows that for the instability do- 
mains either of the two following conditions is satisfied 
in the vicinity of stability-instability boundaries: 
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FIG. 1. Parametric interaction of counterpropagating 
waves, (a) Nonlinear grating is used to couple forward- (" +" ) 
and backward-propagating ("-") waves; (b) square- wave func- 
tion d(z) approximates the periodic modulation of the effec- 
tive second order nonlinear susceptibility coefficient. 

At a+ = a- and S — aa + /2 the system (Eh has an ex- 
act analytical stationary localised solution in the form 



of E±(x) = 3a + /[4|d,|cosh 2 (, 
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For other values of the parameters, analyti- 



cal expressions cannot be found and numerics should be 
used. For numerical analysis it is more convenient to 
renormalise the system (pi) reducing the number of pa- 
rameters. As a result we obtain the system 
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1 + a. The connec- 
1 is given by scaling 
"di\, x = A/ A /oT 



Et = a+V ± I \-l\, -2^--r— +" /l«"«l> ----/ V "+' 
z = Z/a + . System (||) has only two parameters, de- 
fined as 7 = a-/a + , a = 8/a + . Stationary solitons of 
the system (||) can be found numerically, e.g., using the 
relaxation technique, for all 7 > 0, (i w > 0. Some exam- 
ples of stationary solitons of the system rt|) are shown in 
Fig. I 

Criteria (|3|), (EJ) can be rederived for the system 
(||), and then used to calculate the boundary of sta- 
bility area in the (a, 7) plane. For the precise calcula- 
tions of the invariants we use the tangential transfor- 
mation of transverse soliton coordinate X — tan7rA, 



Ref. |T(|. Using this map we cover the infinite inter- 
val in X (—00 < X < 00) by a finite one in X 
(-1/2 < X < 1/2). 
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FIG. 2. Examples of stationary solitons due to counter- 
propagating QPM configuration. Profiles are calculated for 
the system (H) . Solid curves correspond to the forward set of 
waves, dashed lines - to the backward-propagating ones. 
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FIG. 3. (a) Existence and stability diagram for the station- 
ary solutions of the system (a). Solid curves are defined by the 
stability threshold condition (B|) ; (b) , (c) enlarged plots of the 
unstable regions I and II. Soliton profiles calculated for the 
points A and B are shown in Fig. ti. Filled circle corresponds 
to the exact solution (see text). 

To confirm the validity of the stability results given by 
Eqs. (pi) and (0) we analyse numerically an eigenvalue- 
eigenvector problem corresponding to Eqs. (ph linearised 
about stationary solitons of interest. In all analysed cases 
the theoretically predicted stability-instability properties 
were confirmed numerically. However, we should note 
that the theoretical approach which we use does not de- 
scribe so-called oscillatory instabilities (see, e.g., [ pT| ) 
and, although we have not detected such instabilities nu- 
merically, further analysis is necessary to completely rule 
out their possibility. The stability/instability domains 
given by criteria (Q) and (||) are shown in Fig. 3. Direct 
numerical modelling of Eqs. (J5J) confirms the results of 
our stability analysis. Two examples of propagation of 
unstable solitons are presented in Fig. [|. 

Optical power required for generation of counterprop- 
agating solitons of a given beam width can be estimated 
by the method proposed in Ref. |13|. After defining the 
soliton width R s as the maximum width at the half- 
maximum of the the second harmonic amplitude E^, the 
conditional minimum Q m i. n of the total power functional 



Q(a,-y) = R 3 B /2 f+™{\V+\ 2 + \V-\ + 2a\W+\ 2 + 2a\W-\ 2 }dX 
has to be found. Then minimal power density I m i n can 
be calculated as: 
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where nonlinear coefficient Xbulk = ( w2 /c 2 )(2/ceon) 1 / 2 x < - 2 - ) 
is expressed in W~ 1 ' 2 cm~ 1 , x^ is the effective ele- 
ment of the second order susceptibility tensor, L e ff is 
the width of a waveguide. For the first order QPM, 
i.e. I = 1, and the values from Ref. [J5J: tq = 20/im, 
A = 2ttc/uj ~ 1.064/im, refractive index n ~ 1.79, 
2 ) « 6p?7j/R and a waveguide with effective width L e // 
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high. However, experimental progress in quantum well 
technology (see, e.g., 14 1) may make lower order coun- 
terpropagating QPM experimentally possible in the near 
future. 

In conclusion, we have demonstrated the existence of 
solitons that are due to counterpropagating QPM in 
quadratic media. We obtained an analytic criterion for 
the stability threshold for these solitons and found a sub- 
stantial region of stability with only two small regions of 
unstable solitons. We also discussed the conditions for 
experimental observation of these novel solitons. 
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FIG. 4. Examples of unstable propagation of a slightly 
perturbed stationary soliton shown in Fig. 2, A. Two typi- 
cal kinds of instability correspond to slightly increased, (a), 
and reduced, (b), soliton amplitudes. Unstable soliton with a 
slightly increased amplitude is clearly evolving to some stable 
state with an excited internal mode (see, e.g., ]12|). 

Our analysis shows that the point of optimal gener- 
ation is in the domain of stability and corresponds to 
V + = V~ and W + — W~, i.e., one can generate the 
whole four-wave soliton using only two seeded forward- 
propagating waves at one end of a crystal and a mirror 
on the other [the use of mirror will also halve the gen- 
eration threshold (g)]. Thus the solitons due to counter- 
propagating QPM, in principle, require less optical power 
for an experimental observation in comparison with con- 
ventional quadratic solitons, for which the corresponding 
value is I m in ~ 3.4 W/fim (see Ref. [jl3[). In practice 
lowering the generation threshold requires a very short 
(sub-micron) grating period to arrange for lower order 
QPM. This was not the case for the experiments [pj where 
the grating period was about 3/im and QPM order was 
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